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£3 . Exact distribution of the moment estimator of shape parameter for the 

gamma distribution for small samples is derived. Order preserving properties 
' of this estimator are presented. 
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1 Introduction and preliminaries 

It is well known that the gamma distribution has wide application, for example in 
J> \ meteorology to describe the distribution of rainfall. 

qq \ The moment estimators a and A of the parameters a and A of the gamma 

distribution with density 

od ■ i , 

O" fix; a, A) = . - x a exp(— x/X), x > 0, a>0,A>0 

r(a)A Q 



X 



for a random sample X\ , . . . , X n are 



X 2 » S 2 

(i.i) & = ~s 2 and x= x'' 

where X and S 2 are sample mean and sample biased variance respectively. 

In this paper we are mainly interested in order preserving property of the 
estimator a. In general, we say that the estimator 9 based on the sample X 1; . . . , X n 
from population with density fi-;0) is increasing in with respect to the order -< 
if 9i -< §2, whenever 9\ < 9 2 , where 9\ (§2) is the estimator based on sample from 
density /(•; 9{) (/(•; 9 2 ))- We shall consider when -< is one of the following orders: 
<st, <disp, <*• For details and definitions of these stochastic orders we refer the 
reader to Shaked and Shanthikumar [S]. 
Nowak [3] proved the following theorem. 
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Theorem 1. The estimators a and A are stochastically increasing respectively in 
a and A. 



Deriving exact distributions of moment estimators a and A is very tedious for 
n > 2. Now we find the distribution of at for n = 2. First we prove the more general 
fact. 

Lemma 1. Suppose that X\ ~ G(A,ai) and X 2 ~ G(X,a 2 ), X\ and X 2 are 



Corollary 1. Let X\,X 2 ~ G(X,a), Y ~ Ex(l), X\, X 2 and Z are independent 
random variables. Then 



Problem 1. Equation U.fy) yields some interesting property of the exponential dis- 
tribution. There arises a question: Does this equation characterize the exponential 
distribution? 

Corollary 2. For n = 2 the statistics 1/a has the beta distribution B(l/2,a). 
Hence a has monotone likelihood ratio. 

Now we examine dispersive ordering of the a estimator. First we notice that 
1/a is not dispersive monotone in a. We conclude it from Theorem 3.B.14 of 
Shaked and Shanthikumar [5] since the support of the estimator 1/a is the finite 
interval (0, n — 1). 

Theorem 2. The estimator 1/a is not dispersively monotone in a. 

We now consider dispersive ordering of the estimator a. From previous theorem 
it does not follow that a is not dispersive monotone in a, though the function 1/x 
is decreasing and convex for x > 0. We should notice that Theorem 3. B. 10(b) 
of Shaked and Shanthikumar |5J is not valid here. For example, in the following 
theorem we prove that for n = 2 the estimator a is dispersively increasing in a. 

Theorem 3. For n = 2 the estimator a is dispersively increasing in a. 

Proof. Let f a denotes density of the estimator a. We know that 1/a has the beta 
distribution B(l/2,a), so 




(1.2) 



Xf + Xj 

(x, + x 2 y 



(X 1 + Y) 2 + Xj 

(x 1 + r + x 2 y 



i 



3/2 (l - l/x) 



5(1/2, a) 



■x 



X > 1. 
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In order to prove that &\ <disp &2 whenever a,\ < a 2 it suffices to prove that 

S'ifc'ix + 1) - f a2 (x + 1)) < 2 for every c> 

with the sign sequence being — ,+,— in the case of equality and 6t\ < st a 2 (see 
Theorem 2.6 and Remark 2.1 of Shaked [4J). 
Fix c > 0. When x > c then 

f?{ X + 1) - r(x + 1) = 1 (x + i)- 3/2 (i - i/(x + i))-*-^ 

x + l) 3/2 (a; - c + l)- 3/2 (l - l/(o; - c + l)) ai_1 /(l - l/(x + l))" 2 " 1 - l] , 

where A = B(l/2,a 2 )/B(l/2,a 1 ). Thus S^/" 1 - /° 2 ) = S-(h{x)), where 

= log A + 3/2 log(x + 1) - 3/2 log(x - c + 1) + 

+ («! - 1) log(l - l/{x - c + 1)) - {a 2 - 1) log(l - l/(x + 1)). 

Now we show that /i has exactly one maximum when a,\ > 1 and the sign sequence 
is —,+,—. In the other case h is decreasing. After calculating we have that h'(x) = 
if and only if — w(x)/(2x(x + l)(x — c)(x — c + 1)) = 0, where 

= (2a 2 -2ai + 3c)x 2 + (2«2-2ai+4c-4a2C-3c 2 )x + 2c-2a 2 c-2c 2 + 2a2C 2 . 

So we must show that w has only one root in the interval (c, oo) if a± > 1 and no 
roots when «i G (0,1]. Define u)(x) = w(x + c). Therefore 

w(x) = (2(a 2 - a 2 ) + 3c)a; 2 + (2(a 2 - ai) + 4c(l - a x ) + 3c 2 )a; + 2(1 - a x ){c + c 2 ). 

It is easy to see that for ol\ > 1, w has two different roots x\ and x 2 that a;ia; 2 < 
and for ol\ G (0, 1] w has no roots in (0, oo). At the end we must show that the 
sign sequence is — , +, — when S~(h) = 2. For a± > 1 

lim h(x) = — oo, lim h(x) = log A < 0. 

Combining these above facts with Theorem [T] we end the proof. 

Remark 1. The above theorem does not hold in general for n > 3, but direct 
calculating is impossible due to occurrence of hyper elliptic integrals. For example, 
numerical calculations show, that for n = 3, a± ^disp «2 if ai = 1/5 and a 2 = 1/4. 
Then G~ x (x) —F~ l (x) has the local maximum at x ~ 0.72 and the local minimum 
at x w 0.85. 

Deriving exact distribution when n = 3 is more tedious. We should add, that 
the estimator a is closely related to the student statistic y/nX/S. The properties 
of this statistic was very intensively studied in the literature, especially under non- 
normal conditions. First, we find joint distribution of the vector (X, S). One can 
proof the following lemma, see for example Craig [2]. 
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Lemma 2. If f is a density on (0, oo), then the joint distribution (X, S) forn = 3 
is given by formula 



F(x, s) 



18s y(*i)f( 3J ^f ±E )f( 3J ^)dxu 0<s< xV2/2, 

18s 



3x-^6s 2 -3x 2 - , [7L 

f 2 _|_ r^+«v2 

Jo ' Ji 



3x+y , 6s 2 - , .ix 2 
2 

l/( a:i )/(3£z|l±fl)/(Mz| l =fl 



xv/2/2 < s < xv^, 



where R = a/6s 2 — 3(rri — re) 2 . 



If we have the joint distribution (X, 5) it is easy to derive the distribution of 
the statistics T = X / S = \f& from the formula 

fr(t) — J uF(ut,u)du. 

Proposition 1. Forn = 3 the cumulative distribution function ofT for the expo- 
nential distribution with density f(x) = e~ x , x > is given by 



1 - 



F(t) = 



2tt 
3%/3t 2 



V2-t 2 _| 7T_ 



2 arcsin 



On the other hand it is not easy to derive in the general the distribution of T. 
For example, if a G (0, 1) above integrals can be calculated only by numerical 
methods. For n — 3, after a bit algebra we can also derive the distribution of T 
for the gamma distribution with shape parameter a — 2. 

Proposition 2. Forn = 3 the cumulative distribution function of T for the gamma 
distribution with density g(x) = xe~ x , x > is given by 



G(t) = 



10tt(4^-3J 
27v / 3* 4 



t>y/2, 



V2-t 2 (-33t 4 -13t 2 +8)+57rt(4t 2 -3)-30t(4t 2 -3) arcsin (^=) ^ 



27\/3t 6 1 2 

From previous propositions it follows, that the statistic T is not monotone 
with respect to the star order <*, that is the function G~ 1 F(x)/x is not monotone 
in x > 0. To see it, let us calculate: 



x 



1) 

2 ^ G^Fjx) 
3 ^ G-^Fjx) 



10 

33 



\ 



11 ^220tt+ y IOtt (- 891 a/79- 560tt+16200 arccos ( 7^75 ) ) j 
33 v / 79+2007t-600 arccos^Y^) 



32686 



-VI 



=v / 2 



40^+2^10^(2271-27^) 



27^+1871 



1.31502 



V / 40+2v / lf5 
6 



1.32081 
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Since the star order is preserved under increasing function, we have the follow- 
ing corollary. 

Corollary 3. For n = 3 the estimator a is not increasing in a with respect to the 
star order. 

The star order is closely related to the dispersive order, since 

X<*Y iff logX< disp logF. 
Using the same technique as in Theorem 3 we can proof the following proposition. 

Proposition 3. For n = 2 the estimator at is increasing in a with respect to the 
star order. 

Applying the property that for nonnegative random variables such X < st Y and 
X <* Y implies the ordering X <di sp Y (see, for example, Ahmed et al. [I]) we 
can also deduce from Theorem 1 and Proposition 3 that a is dispersively monotone 
for n = 2. 

References 

[1] A. N Ahmed, A. Alzaid, J. Bartoszewicz, S. C. Kochar, Dispersive and super- 
additive ordering, Advances in Applied Probability 18, (1986), 1019-1022. 

[2] Allen T. Craig, The simultaneous distribution of mean and standard deviation 
in small samples, The Annals of Mathematical statistics, Vol 3, No. 2, 1932, 
p.126-140. 

[3] P. Nowak, Stochastic comparisons of moment estimators of gamma distribu- 
tion parameters, to appear in Applicationes Mathematicae. 

[4] M. Shaked, Dispersive ordering of distributions, J. Appl. Prob. 19 (1982), 
310-320. 

[5] M. Shaked, J. G. Shanthikumar, Stochastic Orders, Springer Verlag, New 
York, 2007. 



5 



